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We report the specific heat and entropy of single crystals of the spin ice compound Dy2Ti2 07 
at temperatures down to 0.35 K. We apply magnetic fields along the four characteristic directions: 
[100], [110], [HI] and [112]. Because of Ising anisotropy, we observe anisotropic release of the 
residual zero-point entropy, attributable to the difference in frustration dimensionality. In the high 
magnetic field along these four directions, the residual entropy is almost fully released and the 
activation entropy reaches Rki2. However, in the intermediate field region, the entropy in fields 
along the [111] direction is different from those for the other three field directions. For the [111] 
direction the frustration structure changes from that of three-dimensional(3D) pyrochlore to that 
of two-dimensional(2D) Kagome-like lattice with constraint due to the ice rule, leading to different 
values of zero-point entropy. 

PACS numbers: 65.40.Ba, 65.40.Gr, 75.25,+z, 75.50.Lk 



Recently, there has been increasing attention in the 
physics of geometrical frustration in pyrochlore oxides 
A2B2O7. The spin ice behavior has been observed in 
Ho 2 Ti 2 7 IE II H El, Dy 2 Ti 2 7 0, 1 1 EO, EH and 
Ho 2 Sn 2 07 [13 • Since there is no long-range magnetic 
ordering of the rare-earth moments in these materials at 
least down to 50 mK 0, , the ground state is believed 
to be macroscopicall y d egenerate. In fact, the observed 
residual entropy TT^ shows an excellent agreement 
with the expected zero-point entropy of (l/2)i?ln(3/2) 



In these materials, the A-site ions constitute a three- 
dimensional (3D) network of corner-shared tetrahedra 
(the pyrochlore lattice). Because of the crystal field ef- 
fect, the magnetic moments of the A-site ions, such as 
Dy 3+ and Ho 3+ , have Ising anisotropy along the local 
(111) direction, which points to the center of the tetrahe- 
dron from a vert ice. Owing to this Ising anisotropy, the 
spin responses to magnetic fields are very anisotropic. 
In measurements on polycrystalline samples, there are 
specific-heat peaks at field-independent temperatures of 
0.34 K, 0.47 K and 1.12 K; it was speculated that these 
peaks may be due to ordering of spins with their Ising 
axes perpendicular to the field Q. Moreover, owing to 
the difference in the spin configurations in fields along 
different directions, the process of releasing the resid- 
ual zero-point entropy should be qualitatively different 
reflecting the structure of frustration. However, up to 
present there are few studies of the entropy release using 
single crystals |l3|, Ujj, Hfl • 

In this paper, we report the specific heat and entropy of 
single crystalline Dy 2 Ti 2 07 in magnetic fields along four 
characteristic directions, [100], [110], [111] and [112]. To 
the best of our knowledge, relevant data for the [110] and 
[112] directions have not been published previously. We 
indeed observed an anisotropic process of releasing resid- 
ual entropy for fields along different directions. We will 
discuss this anisotropic behavior in terms of a qualitative 



difference in the dimensionality of frustration structures. 
The single crystals of Dy 2 Ti 2 07 and Eu 2 Ti 2 07 used in 
this work were grown by a floating zone method 8]. We 
measured the specific heat between 0.35 and 40 K and 
in fields up to 5 T by a relaxation method using a com- 
mercial calorimeter (Quantum Design, model PPMS). To 
obtain reliable entropy values, it was essential to measure 
the specific heat of the sample holder (addenda) at each 
corresponding magnetic field. In this way, our data pre- 
sented here does not involve data manipulation to com- 
pensate for the field dependence of the background. 
The sample size was approximately 2.0 x 2.0 x 0.06mm 3 
and the mass was about 1.7 mg. In order to avoid field 
inhomogeneity due to the demagnetization effect, we cut 
the crystal along the (110) plane. All the four character- 
istic directions of the cubic lattice, [100] (equivalent to 
[001]), [110], [111] and [112] (Fig. ffl (a)), lie in the sur- 
face plane of the plate-like sample. This configuration 
allows us to examine the anisotropy based on the data 
taken from the same crystal. 

Because Dy 2 Ti 2 07 is an insulator, the total specific heat, 
C-totai? can be expressed as 
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where C sp i n and C pn0 non represent electronic-spin and 
phonon contributions, respectively. The nuclear Schot- 
tky term is negligible in this compound. In order to esti- 
mate the phonon contribution, we also measured the spe- 
cific heat of Eu 2 Ti 2 07. Because of the similarities of the 
lattice parameters and the atomic mass of the A-site ions, 
and of the nonmagnetic ground state and insulating prop- 
erty of Eu 2 Ti 2 07, we approximate C pn onon of Dy 2 Ti 2 C>7 
with the total specific heat of Eu 2 Ti 2 07 with a small 
correction for the difference in A-site mass, {M^ y /M^ u ) 
C t / otal (Eu 2 Ti 2 07). Thus the spin contribution is ap- 
proximated as C sp in(Dy 2 Ti 2 C>7) = C t otai(Dy 2 Ti 2 7 ) - 
(M Dy /M E u)C t / otal (Eu 2 Ti 2 7 ). Here M Dy and M Eu rep- 
resent molar mass of Dy 2 Ti 2 07 and Eu 2 Ti 2 07, respec- 
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tively. We estimate the activation entropy S a (T) by in- 
tegrating C spin /T from 0.35 K to T of up to 40 K: 

S a (T) = f %^dT (2) 

J 0.35 K 1 

We deduce the residual entropy, So , as Rlu.2 — S a (T\ ) , us- 
ing the activation entropy at high temperature T\ ^> J e ff 
where J e ff is the effective interaction between the nearest- 
neighbor spins. In Dy2Ti2 07, the main contribution to 
Jeff = 1-11 K is the dipolar interaction and is ferromag- 
netic (fm) nmos. 

In Fig. [21 (a) we show the spin component of the specific 
heat measured after field-cooling at 0.5 T along the four 
characteristic directions. It is estimated that the first ex- 
cited level under the crystal field is nearly 150 K [HI, so 
that the total spin entropy involved at low temperature 
is RYn2. The Zeeman energy of the applied field for flip- 
ping a spin is expressed as AEz = 2gjJ hb{^oH) x cos#, 
where gj is the Lande g factor, J is the total angular mo- 
ment of Dy 3+ ion 19], /lib is the Bohr magneton, and 
is the relative angle between the local (111) and the field 
direction. The Zeeman energy for a field of n$H = 0.5 T 
along the Ising axis is 6.8 K. In a pyrochlore lattice, the 
four spins in each tetrahedron have Ising axes along lo- 
cal (111) directions and thus point to different directions; 
hence the average value of the Zeeman energy is smaller 
than this value. In fact at 0.5 T, the Zeeman energy for 
one of the four spins exceeds the internal field governing 
the ice rule, 4J e ff. It should be noted that Fig. (b) 
clearly shows that only for the [111] field direction there 
still remains residual entropy of So « 0.44 J/K mol-Dy. 
Let us first examine the origins of the specific heat peaks 
based on the spin ice model [l3|. In fields along the [100] 
direction, the stable spin configuration is uniquely de- 
termined as two-spins in and two-spins out (2-in 2-out), 
satisfying the ice rule. Moreover, all the spin compo- 
nents parallel to the field are equivalent. Therefore there 
should be only one peak in C sp i n (T)/T associated with 
the thermal excitation across the energy gap determined 
by 1 4 J e ff + AEz\. Since cos# is 1/y/S the additional Zee- 
man energy term at 0.5 T is 3.9 K. This quantitatively 
explains the observed linear increase of the peak temper- 
ature with the field strength in a high field region [l3j . 
In fields lower than 2 T, a changeover to behavior with a 
nearly field-independent specific heat peak attributable 
to J e ff is observed. 

In the [110] field direction, two of the four spins in 
each tetrahedron is determined as 1-in 1-out with cos# 
= <s/2/<\/3, giving the expected additional Zeeman en- 
ergy term of 5.4 K. The other two spins are perpendic- 
ular to the applied field; their directions are not deter- 
mined by the external field. Thus the former spins should 
give rise to a Schottky peak with a peak temperature 
increasing with the field, whereas the latter spins may 
lead to specific heat anomaly associated with ordering 
of these sublattice spins due to their interactions. Since 
the peak position shown in Fig. [3 (a) is unchanged with 
the field strength for H\\ [110], we attribute the field- 



independent peak at 1.1(1) K observed in polycrystalline 
samples 0, [l3| to the ordering of these two free spins. 
In the [111] field direction, one of the four spins is parallel 
to the field and most effectively stabilized by the exter- 
nal field. However for the other three spins, the nearest 
neighbor interaction dictating the ice rule competes with 
the Zeeman energy of the external field, because the lat- 
ter demands the 1-in 3-out configuration. Therefore the 
energy to flip one of these three spins is the difference be- 
tween those two interactions. The energy gap becomes 
1 4 J e ff - AEz\ with cos# = 1/3, which reaches zero at 1.0 
T. In lower fields, the stable spin configuration is the 
ice-rule 2-in 2-out state because it is dominated by the 
spin- spin interaction. However, in higher fields the ice 
rule breaks down and the stable configuration changes 
to 1-in 3-out. This crossover process has been quantita- 
tively confirmed in the magnetization measurements 
In another characteristic field direction, the [112], the 
configuration of three of the four spins are specified as 1- 
in 2-out. The value of cos# is 2\j2jZ for one and y/2/3 for 
the other two, giving the expected additional Zeeman en- 
ergy at 0.5 T of 6.4 K and 3.2 K, respectively. The other 
spin in a tetrahedron is perpendicular to the field. Be- 
cause of the internal interaction, however, this spin points 
outward to satisfy the ice rule. Therefore the energy dif- 
ference between up and down states for this spin depends 
not on the external field strength, but only on the internal 
field. There is a possibility that this field-perpendicular 
spin is the origin of some of the other field-independent 
peaks observed in polycrystalline specific heat. However, 
for the [112] field direction we observed only one peak, 
which depends on H (data not shown). 
Let us now compare in Fig. (b) the spin contribution 
to the entropy S a at 0.5 T along the four characteristic 
directions and that of a single crystal at T. Let us first 
examine the observed entropy release at low tempera- 
tures. In the low temperature region, the entropy is most 
actively released for the [111] field direction, because the 
energy gap is reduced by the external field as explained 
above. In the other characteristic field directions, exter- 
nal fields do not compete with the ice rule; they stabilize 
some particular spin configurations within the spin-ice 
manifold. It is thus the number of field-perpendicular 
spins which mainly determine the activation entropy at 
low temperatures, since these spins are not further sta- 
bilized by the magnetic field. The numbers of such spins 
are two, one, and zero for H // [110], [112], and [100], 
respectively. The activation entropy at low temperatures 
indeed decreases in this order as shown in Fig. 0(b). 
Next, we focus on the values of the activation entropy in 
the high temperature region. The difference from Rhi2 
represents the residual entropy, which is not released by 
the external field, and are intimately related to the frus- 
tration structure. In the [111] field direction, the near- 
est neighbor interaction governing the ice rule is stronger 
than the external field for three of the four spins up to the 
field of 1.0 T. Thus at 0.5 T the frustration still remains. 
These three frustrated spins constitute a Kagome lattice 
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(Fig. 3). Viewed from the [111] direction, the pyrochlore 
lattice consists of a layered stacking of a triangular lattice 
and a Kagome lattice (a pyrochlore slab). The frustra- 
tion on a 3D pyrochlore changes to that on a 2D Kagome 
lattice. The ground state of an ordinary Kagome lattice 
has grea ter residual entropy than that of pyrochlore lat- 
tice [2(1 |^J. However, in the present case the residual 
entropy of this Kagome lattice is smaller than that of a 
pyrochlore lattice because spins on this Kagome lattice 
have the Kagome-ice rule constraint that the "up trian- 
gle" satisfies 2-in 1-out and the " down triangle" satisfies 
1-in 2-out (Fi g. 3) . Thus, this state may be called as 
"Kagome ice" [22j. Its residual entropy is estimated as 
|jRln3ir(!)£ =. 0.598 J/K mol-Dy by Pauling's method 
[23q . Udagawa et al. calculated the exact value for the 
Kagome ice model to be 0.672 J/K mol-Dy |24|. 
As shown in Fig. 2(b), in fields along the other three di- 
rection the residual entropy is fully released (5 a = 5.78 
± 0.08 J/K mol-Dy) but along [111] direction, S a = 5.35 
J/K mol-Dy. Thus in this study the residual entropy is 
evaluated as So = 0.44 ± 0.12 J/K mol-Dy. This value 
is somewhat smaller than the theoretical expectations. 
At 1 T along the [111] direction, the effective Zeeman 
energy overcomes the energy scale of the nearest neigh- 
bor FM interaction and spin configuration becomes the 
ice-rule breaking 1-in 3-out. Therefore the residual en- 
tropy should vanish and this expectation agrees with our 
experimental results (Fig. 4). Note that in refs. [l5lH^. 
they reported for the Kagome-ice residual entropy of 0.8 
J/K mol-Dy and 0.65 J/K mol-Dy after subtracting the 
extrinsic background estimated by a quadratic fitting. 
In zero field, the role of long-range dipolar interaction is 
known to be quite important for quantitative description 
of the spin ice ^tJ. F° r the Kagome ice realized in the 
[111] field direction, its lower dimensionality should make 
less important the long-range dipolar interaction within 
the Kagome lattice, since the increase in the number of 
lattice sites with the distance is not as rapid as to com- 
pensate for the decay in the dipolar interaction. We note, 
however, that the dipolar field from the ordered moment 
of the triangular lattice provides an additional molecular 
field not included in the " Kagome-ice rule" . Such dipo- 
lar field at a Lattice-lattice site amounts to a correction 
of -19% to the nearest-neighbor contribution if the con- 
tribution from up to the fourth-nearest-neighbor sites on 
the triangular lattice is included. 

In the characteristic field directions other than [111], the 
zero-point entropy is more readily released. In the [100] 
and [112] field directions, the spins are uniquely deter- 
mined in the 2-in 2-out configuration without degener- 
acy. In the [110] field direction, the structure of frus- 
tration changes to a one-dimensional (ID) chain without 
macroscopic degeneracy as shown in Fig.l (b). In real- 



ity, to overcome the spin freezing setting in below ~1K 
external field of ~ 0.1 T is needed to remove the spin-ice 
residual entropy even in the field-cooling process. Fig. 2 
(b) indicates that 0.5 T is indeed sufficient to remove the 
residual entropy. 




FIG. 1: (a) Four characteristic directions in which the mag- 
netic field is applied, (b) ID chain of spins perpendicular to 
the field on a pyrochlore lattice in the [110] field direction. 



In conclusion, we presented anisotropic release of the 
zero-point entropy of single crystalline Dy2Ti2 0y, as- 
cribable to the qualitative difference in the structure of 
frustration. In the [111] field direction, we consider this 
system as a dipolar Kagome ice, in which the spins frus- 
trate on a 2D Kagome lattice with ice-rule constraint 
|25j in the presence of additional dipolar field from the 
ordered moments on the triangular lattice. We estimate 
the residual entropy by considering the frustration only 
on the Kagome lattice, but the role of the dipolar field 
from the triangular-lattice slab is expected to be impor- 
tant in the value of the residual entropy, as well as in 
the process of entropy release. Such difference between 
the Kagome-ice behavior and the "dipolar Kagome-ice" 
behavior is worth further attention in the future studies. 
We observed a nearly field-independent peak in the spe- 
cific heat in fields along the [110] direction. This is most 
probably the origin of one of the field-independent peaks 
observed in polycrystalline samples at 1.1 K 0, In 
the [112] field direction, one of the four spins is perpendic- 
ular to the external field. Although there is a possibility 
that ordering of this "free" spin is the origin of another 
field-independent peak observed in polycrystals, we did 
not observe any sign of such a peak in the present mea- 
surements down to 0.35 K. 
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FIG. 2: The spin contributions to (a) specific heat and (b) ac- 
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and negative, respectively, by the ice rule. 
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